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' Abstract. We report a theoretical study of sulphur chemistry, as applied to hot 

cores, where S-bearing molecular ratios have been previously proposed and used 
as chemical clocks. As in previous models, we follow the S-bearing molecular com- 



position after the injection of grain mantle components into the gas phase. For 
this study, we developed a time-dependent chemical model with up-to-date re- 
action rate coefficients. We ran several cases, using different realistic chemical 



compositions for the grain mantles and for the gas prior to mantle evaporation. 
The modeling shows that S-bearing molecular ratios depend very critically on 
the gas temperature and density, the abundance of atomic oxygen, and, most im- 
portantly, on the form of sulphur injected in the gas phase, which is very poorly 
known. Consequently, ratios of S-bearing molecules cannot be easily used as chem- 
ical clocks. However, detailed observations and careful modeling of both physical 
and chemical structure can give hints on the source age and constrain the mantle 
composition (i.e. the form of sulphur in cold molecular clouds) and, thus, help 
to solve the mystery of the sulphur depletion. We analyse in detail the cases of 
Orion and IRAS16293-2422. The comparison of the available observations with 
our model suggests that the majority of sulphur released from the mantles is 
mainly in, or soon converted into, atomic form. 
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1. Introduction 



It is a long-standing dream to use relative abundances of different molecules as chemical 
clocks to measure the ages of astronomical objects. Studies of t he ag es of star formation 
region s have recently focused on S-bearing molecules 



Charnlev 




and 



Hatchell et al 



1 1998) were the first to propose that the relative abundance ratios of SO, SO2 and H2S 
could be used to estimate the age of the hot cores of massive protostars. The underlying 
idea is that the main reservoir of sulphur is H 2 S on grain mantles, and that when the hot 
core forms, the mantles evaporate, injecting the hydrogen sulphide into the gas phase. 
Endothermic reactions in the hot gas convert H2S into atomic sulphur and SO from which 
more SO and, subsequently, SO2 are formed, making the SO2/SO and SO/H2S ratios 
nice functions o f time. These studies have triggered a variety of work, both observational 



and theoretical J Charnlev 



1997 



Hatchell et al 



1998: 



Buckle fc Fuller 



2003). 



This line of research, however, has been challenged by recent ISO observations, which 
have cast doubt on the basic assumption that sulphur is mainly trappe d in grain mantles 
as H 2 S. The lack of an approp riate feature in the ISO spectra of high IjGibb et al.ll2000n 



and low IjBoogert et al 



2000) mass protostars sets an upper li mit on the ma ntle H2S 



abun dance which cannot exceed about 10~ 7 with respect to H2 Ijvan Dishoeck fc B lake 
1998|) . Indeed, the identity of the major reservoir of sulphur in cold molecular clouds 
is a long standing and unresolved problem, for the sum of the dete ctable S-bearing 



molec ules is only a very small fraction of the elemental S abundance (Tieftrunk et al 



1994). Since sulphur is known not to be depleted in the diffuse medium fe.g. ISofia et al 



19941) , it is usually assumed that sulphur in dense clouds is depleted onto the grain 



mantles rather than in refractory cores (e.g. ICasell 



Hasegawa 



Ruffle et al 



& H erbst 



1994), but 



1999) proposed that 



how this happens is a mystery. In a theoretical study, 
in collapsing translucent clouds sulphur is efficiently adsorbed onto grain mantles. In 
fact, in these regions, most of the gas-phase sulphur is in the form of S + , while grains are 
typically negatively charged, so that the collisional cross section for sulphur is enhanced 
compared with neutral species (e.g. O) and sulphur is removed from the gas phase more 
rapidly. 

Another mystery is the form of sulphur on dust grains. The simplest possibility is 
that it consists of relatively isolated atoms, as would occur in a matrix, or perhaps as 
isolated pairs of atoms (S2). Another possibility is that the sulphur is amorphous (or 
even crystalline), having formed islands of material from the initially adsorbed atoms. 
Crystalline sulphur is known to come in two forms - rhombic and monoclinic - both of 
which consist of Ss cyclic molecules. Vaporization leads to a complex mixture of sulfur 
polymers through Ss in complexity. If sulphur is elemental and amorphous, evaporation 
is also likely to lead to molecules of sulphur through eight atoms in complexity. So far, the 
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only S-bearing species firmly detect ed on granular surfac es is OCS, but with a relatively 

1 20021) 



low fractional abundance of 10 7 (Palumbo et al 



1997). Recently, 



Keller et al 



claimed the detection of iron sulphide (FeS) grains in protoplanetary disks, but there is 
no evidence to suggest that solid FeS is the main form of sulphur in the parent collapsing 
environment. Actually, if th e main form of so lid sulphur is FeS, S s hould follow Fe deple- 



tion, which is not observed iSofia et alJl!994fl . Even more recently, 



Scappini et al 



(2003) 



suggested that hydrated sulphuric acid (H2SO4 • H2O) is the main sulphur reservoir. 
In whatever form sulphur resides in the grain mantles, there is the possibility that the 
species, once evaporated, are very quickly destroyed to give atomic sulphur. In summary, 
although all the evidence is that sulphur is depleted onto grain mantles in cold clouds, 
its particular form is very uncertain. 

Given the need for chemical clock methods, it is timely to reconsider the use of 
S-bearing molecules in this fashion. In this paper, we present a model with an up-to- 
date chemical network involving S-bearing molecules. We run several cases to cover a 
large, realistic parameter space for hot core sources, consistent with present observational 
constraints. Based on the results we obtain, we conclude that it is tricky to use abundance 
ratios of S-bearing molecules as chemical clocks in the absence of other constraints, for 
they depend more on the initial conditions, gas density, temperature, and the initial form 
of sulphur injected in the gas phase than on the age of the source. 

The paper is organized as follows: we describe the model in §2, the model results in 
§3, and in §4 we discuss the practical consequences of those results and apply the model 
to the specific cases of Orion and IRAS16293-2422. 



2. The Model 

We have developed a pseudo-time dependent model for the gas phase chemistry that 
computes the evolution of the chemical composition of a volume of gas with a fixed den- 
sity and temperature. Our goal is to follow how the S-bearing molecular abundances vary 
with time when the gas undergoes a sudden change in its temperature and density, and/or 
in its overall chemical abundance, because of the evaporation of grain mantles. In hot 
cores the dust temperature increases to an extent that it exceeds the mantle evaporation 
temperature, i.e. ^100 K, and all the components of the grain mantles are suddenly in- 



(e.g. 



Brown et al. 


1988 


Charnlcv et al. 


1992; 


Casclli et al. 


1993 


Millar et al. 


1997b) 



fact, it is more probable that hot cores have spatial gradients in temperature and density 
relea sing the molecules at different times depending on their surface binding energies 
( ^eccarelli et al.lll996 ; 



sec 



Viti fc Williams 1999; 



.1 



Dotv et al.H2002 ; 



Rodgers fe Charnlev 



20031) . However, the goal of this work is mainly to test the effects of the form of the main 



initial sulphur bearing molecules and we preferred to simplify the problem assuming that 
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all the molecules evaporate simultaneously from the grain. More detailed models will be 
presented in a forthcoming paper. 

In order to simulate these conditions, the gas-phase chemical composition prior to 
evaporation of the mantles is taken to be similar to that of dark molecular clouds. At 
time t= 0, the grain mantle components are injected into the gas phase, and the model 
follows the changes in the gas chemical composition with a given gas temperature and 
density. Throughout this paper we will use the word "evaporation" to refer to the loss of 
the grain mantles. 

The model is a reduced chemical network, which includes 930 reac- 
tions involving 77 species containing the elements H, He, C, O and S. 
The standard neutral-neutral and ion-neutral reactions are considered. Most 
of the reaction coefficients are from the NSM ("new standard model"; 
http://w ww.physics.ohio — state . edu/ eric /research- files /cddata.julyOS) database; 



see also 



Lee. Bettens. fc Herbstl l|199fil) . updated with new values or new analy- 



ses of assorted values in databases (e.g. the NIST chemical kinetics database at 



http://kinetics.nist.gov/index.php) when available. Furthermore, several high temper- 



ature (neutral 
have followed 



neutral) reactions have been added. To select the reduced networ k, we 



and 



Ruffleet al 
fl) 



(2002) for CO formation, 



Hollenbach & McKcc (197 93) and 



Hartauis t et al.l (1198011 for the oxygen chemistry and iPineau Pes Forets et al 



(1993) 



Charnlevl l|1997|) for the sulphur chemistry. 



To vali date this network, we c o mpar ed our results with abund ances prev i ously ob- 



tained by 



Lee. Bettens. &: Herbstl 1)1996|) at low temperature and 



Charnlevl l|l997j) at 



higher temperatures using the same initial abundances as these authors. We found that 
we can reproduce molecular abundances to better than a factor of three. This is an 
indication that small variations in the rate coefficients between the updated NSM and 
UMIST databases do not strongly influence the computed abundances of sulphur bearing 
species. One exception is the CS molecule, which we produce at an abundance ten times 
less than Charnley's model, because our adopted rate coefficient for the reaction CO + 
CRPHOT — > O + C, (where CRPHOT is a photon induced by cosmic rays) is 50 times 
smaller than in the UMIST database. The lowered abundance of C then translates into 
a lowered abundance for CS, since C is a precursor of CS. 

Note that we have assumed the gas to be totally shielded from the interstellar UV 
field and no other UV field to be present. Thus, the model does not include any pho- 
tochemistry, with the exception of cosmic ray-induced photodestruction reactions. The 
model takes into account a reduced gas-grain chemistry: H2 is formed on grain surfaces 
and t he recombination of ions with negatively charged grains occurs (see 



Aikaw a et al 



19991 for the recombination of ions with negative grains). Moreover, neutral species 



can deplete onto grain mantles, and mantle molecules can evaporate because of therm al 



effects and cosmic rays l|Hasegawa. Herbst. fc Leung 



1992; 



Hasegawa fc Herbst 



1993). 
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2.1. Gas-phase and mantle abundances prior to evaporation 

To help determine a set of molecular abundances prior to mantle evaporation, we ran 
the model with a temperature equal to 10 K and a density nu 2 equal to 10 4 cm~ 3 , 
including freeze out, for 10 7 yr. At this time, species such as SO, SO2, and CS reach 
abundance s similar to those ob served in dark clouds (<~ 10 -9 , < 10 -9 , and ~ 10 -9 , re- 
Dic kens et al d elemental a bundances (with re spect to H 2 ) 



spectively; 

for He, O, C+ and S+ are respectively: 0.28, 6.38 x 10" 4 i|Mever et alJll998|) . 2.8 x 10~ 4 



1 Cardelli et al 



1996) and 0.3 x 10 -4 ; the sulphur abundance refers to the gaseous and 
grain portions (see below). The late-time abundances obtained could not be reason- 
ably used without modification for the pre-evaporated chemical compo sition for several 
reasons. First, our model, like most other gas-phase treatments (e.g 



Mill ar et al 



Lee et al 



1996 



1997at). overestimates the O2 and H2O abundances in cold dense clouds 



by orders of magnitude with respect to the ISO, SWAS and ODIN observations (e.g. 



Berein et al. 


2003 


Pagani et al. 


2003) 



2003). Secondly, the elemental sulphur abundance per- 
taining to the gas must be lowered to avoid getting very high abundances of sulphur- 
bearing species. The portion of the abundance not used for the gas can be considered 
to reside in grain mantles until evaporation or in grain core. In order to mimic realistic 
conditions, we thus adopted three different compositions for the pre-evaporative gas, as 
follows: 

Composition A: We adopted the computed late-time molecular abundances except for 
O2 and H2O, which were assumed to be 10 -7 and 10 -8 with respect to H2 respectively, 
in agreement with observations in molecular clouds , and for atomic oxygen, O, which 
was assumed to carry the oxygen not locked into CO, leading to a fractional abundance 

Caux et alJl999: 



of 2.6 x 10 \ as suggested by obs ervations (e.g 



Vastel et al 



2002: 



Lis et al 



Baluteau et al. 1997 



200 1|) . In addition, the initial gas-phase sulphur abundance 
was taken to be a factor of 30 lower than the elemental abundance. 

Composition B: We re-computed the late time abundances, lowering artificially by 
two orders of magnitude the rate of the dissociative recombination of H30 + , to decrease 
the computed O2 and H2O abundances. In this case, it was only necessary to lower the 
initial elemental sulphur abundance by a factor of five. The abundance of atomic oxygen 
in this case is 4.8 x 10~ 4 , consistent with observations in molecular clouds. 

Composition C: The abundances we re taken to equal those measured in the direc- 
tion of L134N, as reported in Table 4 of lCharnlev. Rodeers. fc Ehrenfreundl 1 2001). The 



oxygen not contained in the species reported in Table 4 was assumed to be in atomic 
form. 

All three gas-phase compositions have large abundances of atomic oxygen, in agree- 
ment with observations. We assume that this large O abundance is also present at the 
beginning of the hot core phase, in contrast with previous studies. This implies many 
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Table 1. Computed late time (A, B) and adopted (C) gas-phase abundances with respect 
to H2 for the main S-bearing molecules prior to mantle evaporation. 



Species 


A 


B 


c 


SO 


9 x 10~ 


-9 


3 x 10" 


-9 


3.1 x 10" 9 


SO2 


4 x 1(T 


10 


1.5 x 10 


-10 


1 x 10~ 9a 


H 2 S 


1 x 1(T 


10 


1 x 10~ 


11 


8 x lO" 10 


CS 


4 x 10" 


10 


2 x 10" 


-8 


1.7 x 10" 9 



a Here we took the value of the upper limit. 



Table 2. Adopted abundances for mixtures of evaporated S-bearing molecules (with 
respect to H 2 ). 



Species 


Mod. 1 


Mod. 2 


Mod. 3 


Mod. 4 


OCS 


10- 7 


10" 7 


10- 7 


10" 7 


H 2 S 


10- 7 


10" 7 


10- 7 


10" 8 


S 





3 x 10~ 5 





3 x 10~ 5 


S 2 








1.5 x 10 -5 






differences in the computed abundances, as shown in Section 3.1 and discussed in Section 
4.1. Tableniists the abundances of the main S-bearing molecules for the three gas-phase 
compositions prior to mantle evaporation. 

At time t = 0, the grain mantle components are injected into the gas phase. The 
abundances of the major mantle components are relatively well constrained by the ob- 
servations and we t ook the observed ab undances (with res pect to H2) in high mass pro- 
tostars: H 9O: 10~ 4 |Sc 
4 x 1CT 6 (Chjaretal 



mtte et al 



1996), H^C O: 4 x 10~ 6 llKeane et al 



1996) and CH 4 : HT 6 (Bo ogert et al 



2001), CH3OH: 



1998). Note that in order 



to shorten the number of treated species, we n egleted the CO? molec ule which is abun- 



Gerakines et al 



1999), because it is not 



dant in mantles (2 x 10~ 5 with respect to H2, 
a crucial element of the sulphur chemistry. On the contrary, the situation is very uncer- 
tain with respect to the S-bearing mantle molecules, as discussed in the Introduction. 
In order to study the influence of the injected S-bearing abundances on the evolution of 
the chemical composition, we have run models with four types of material evaporating 
from mantles that differ in their major sulphur-bearing species, as reported in Tabled 
In practice, sulphur on the grain mantles can be stored as OCS, H 2 S, or pure sulphur 
in a matrix-like, amorphous, or even crystalline form. Of the four mixtures, the first one 
(used in Model 1) has the bulk of the sulphur in the refractory core of the grain, while 
the other three (used in Models 2-4) have large abundances of elemental sulphur leading 
upon evaporation directly or eventually to either S or S2 in the gas. 
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2.2. "Important" reactions 

As the post-evaporative gas-phase chemistry proceeds, it is important to determine the 
reactions that will influence the formation and destruction of sulphur-bearing species the 
most severely. We will call them "important" reactions in the following discussion. By 
this term, we mean quantitatively those reactions that lead to significant variations of 
the main S-bearing abundances when the relevant reaction rate is changed by a small 
amount (specifically 10%). Although such a determination has not been featured in pa- 
pers on astrochemistry, we thought it worthwhile to introduce a suitable procedure here. 
The aim of this study is twofold: (i) to determine a set of reactions to check carefully 
in laboratory experiments because the computed abundances are particularly sensitive 
to those reactions, and (ii) to ascertain whether different chemical networks will lead to 
different results, and why. To identify these reactions, we first defined a "perturbation" 
in the rate coefficient for each of the 750 gas-phase reactions by multiplying them by a 
factor of 1.1, one at a time. For each single perturbation, we then computed the abun- 
dance "variations" AX by comparing the computed abundances after 10 4 yr with the 
reference abundances calculated with the non-modified set of reaction rates, according 
to the expression AX — - — — -, where X re f is the reference abundance and X the 
abundance obtained with the perturbed rate. 

Fig^Jshows the variations of the abundances of the main S-bearing species H2S, CS, 
and SO divided by the amplitude of the perturbation (AR = 0.1) for a temperature of 
100 K and a density of 10 7 cm~ 3 . The calculation has been performed for composition 
A and Model 2. The abscissa consists of the numbers of reactions in our network from 
86 to 834. Vertical lines are included for those reactions that produce a normalized 
variation larger than 1% of the largest variation, while the actual numbers of reactions 
with a normalized variation greater than 5% are listed. Note that for a variation (0.1) 
equal in size to the amplitude of the perturbation, the line extends upward to unity 
and that a AX/ AR of 0.5 implies that the abundance of the studied molecule change 
by 5% in abundance upon a 10% change in the reaction rate. The sets of numbered 
"important" reactions for the molecules H2S, CS, and SO consist of 28, 24, and 17 
reactions, respectively. We do not show the variations of OCS and SO2 because the OCS 
abundance is only affected significantly (9%) if reaction 757 is perturbed, while SO2 
shows the same behavior as SO so the corresponding figure is the same. 

In Tableau we list the 70 most important reactions for the chemistry of H2S, CS, SO, 
SO2 and OCS for initial composition A and the four models (see previous section) at a 
temperature of 100 K and a density of 10 7 cm~ 3 . The first two columns give the number 
of the reaction and the actual reaction, while the third column gives the reaction rate 
coefficient k in terms of the standard parameters a (cm 3 s _1 ), (3 and 7 (K): 

k = a(T I '300 \fexp{--y/T) (1) 
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Fig. 1. The "variation" AX (see text) of the abundances of SO, H2S, and CS divided by 
the amplitude (AR = 0.1) of the perturbation of the rate coefficients plotted for reactions 
with number between 86 and 834; i.e., all the reactions in the model exept those dealing 
with dust grains. The numbers indicate those reactions that induce a variation of more 
than 5% of the highest variation for each species. 



where T is the gas temperature. In this table, we report only the rate coefficients that 
differ by more than 10% from the UMIST ones. In the last four columns, the symbol 
• indicates that the reaction is important for sulphur chemistry in the model of the 
corresponding column. The list of reactions in Table depends only weakly on the choice 
of initial composition mentioned in section 2.1, because the three initial compositions A, 
B and C give sufficiently similar abundances at 10 4 yr. 

Of course, for very different initial compositions, the list of important reactions could 
be slightly different, especially if the network is different. This is for example the case of 



the ICharnl ev (1997) model, where all reactions involving atomic oxygen are not impor- 
tant, for no initial gaseous atomic oxygen is assumed to be present, whereas ion-molecule 
reactions involving molecular ions such as fFj~, H30 + , and H3CO + , are crucial. 

We found that the relative importance of the reactions in Table depends more on 
the initial mantle composition (Model 1 to 4) than on the pre-evaporation composition 
of the gas (A, B or C). For Model 1, the chemical network is simpler (i.e. fewer important 
reactions) than for Models 2, 3, and 4 because there is no initial S or S2, so that many 
reactions forming CS as well as reactions with S and S + lose importance. Models 2 and 3 
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have a similar list of important reactions but some reactions with SO gain in importance 
for Model 3. Models 2 and 4 are even more similar in their lists of important reactions. A 
few reactions with H 2 S become less important for Model 4 because the initial abundance 
of H2S in this model is ten times less than in Model 2 and only one reaction (number 
550) becomes important. 

Table 3: List of the reactions important for the chemistry of the 
SO, SO2, OCS, H 2 S and CS molecules for a temperature of 100 K 
and a density of 10 7 cm~ 3 with initial composition A. 



Reaction 1 



a 







7 



Mod. 1 Mod. 2 Mod. 3 Mod. 4 



H 2 + CRP — » H + + H + e~ 

H 2 + CRP -► H+ + e~ 

Hc + CRP -► Hc+ + e~ 

CH 4 + CRPHOT -» CH 2 + H 2 

H 2 + CRPHOT — ► OH + H 

S + CRPHOT -► S+ + 0- 

H 2 S + CRPHOT -> S + H 2 

H 2 S + CRPHOT -» H 2 S+ + e" 

SO + CRPHOT -> SO+ + e- 

OCS + CRPHOT -> S + CO 

OCS + CRPHOT -> OCS+ + e" 

C + SO -> CS + O 

C + S0 2 -> CO + SO 

CH 2 + S -> CS + H 2 

O + CS -> CO + S 

2 + S -> SO + O 

H+ + O -> OH+ + H 2 

H+ + H 2 -> H3O+ + H 2 

H+ + CH4O -» CH+ + H 2 + H 2 

H+ + S -> HS+ + H 2 

H+ + SO -> HSO+ + H 2 

H+ + OCS -> HOCS+ + H 2 

Hc+ + H 2 -> OH + He + H+ 

Hc+ + CO -» O + C+ + He 

Hc+ + CH4O -> CH+ + OH + He 

Hc+ + OCS -» CS + 0+ + He 

CH+ + O -> HCO+ + H 2 



2.860e-19 0.00 



6.08c-14 
2.52e-14 
2.49e-14 
1.34c-13 
4.41c-14 
1.30c-14 
1.39e-13 
3.76e-14 
3.50c-ll 



8.40e-09 
3.80e-09 
1.32e-09 
5.50e-10 
1.70e-09 
8.40c-10 
2.05e-10 



0.00 
0.00 
0.00 
0.00 
0.00 
0.00 
0.00 
0.00 
0.00 



-0.50 
-0.50 
-0.50 
-0.50 
-0.50 
-0.50 
0.00 



0.00 



0.00 
0.00 
0.00 
0.00 
0.00 
0.00 
0.00 
0.00 
0.00 



1.94e-ll 0.00 231.00 

2.30e-12 0.00 0.00 

4.50e-09 -0.50 0.00 

1.80e-09 -0.50 0.00 



0.00 
0.00 
0.00 
0.00 
0.00 
0.00 
0.00 



1 Only rate coefficients that differ by more than 10 percent from the UMIST value are included. 
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437 


CH^ + -> 


HOC+ + H 2 










447 


cuf + SO - 


+ HOCS+ + H 2 


4.20e- 


09 


-0.50 


0.00 


459 


CH 4 + S+ — > 


H 3 CS+ + H 


1.40c- 


10 


0.00 


0.00 


504 


+ HS+ — > 


S+ + OH 










505 


+ HS+ — > 


SO+ + H 










506 


+ H 2 S+ — 


HS+ + OH 










507 


+ H 2 S+ — 


SO+ + H 2 










537 


H 2 + HCO 


f —> CO + H 3 0+ 


2.10c- 


09 


-0.50 


0.00 


539 


H 2 + HS+ 


-> S + H 3 0+ 










540 


H 2 + H 2 S+ 


HS + H 3 0+ 


7.00e- 


10 


0.00 


0.00 


550 


H 2 0+ + H 2 S 


H 3 S+ + OH 


7.00e- 


10 


0.00 


0.00 


552 


H 3 0+ + H 2 CO H 3 CO+ + H 2 


2.60e- 


09 


-0.50 


0.00 


553 


H 3 0+ + CH 4 CH5O+ + H 2 










556 


H3O+ + H 2 S 


H 3 S+ + H 2 










579 


HCO+ + OCS -» HOCS+ + CO 


1.50e- 


09 


-0.50 


0.00 


585 


H 2 CO + S+ 


-f H 2 S+ + CO 


1.10c- 


09 


-0.50 


0.00 


586 


H 2 CO + S+ 


-► HCO+ + HS 


1.10c- 


09 


-0.50 


0.00 


587 


H 2 CO + H 3 S 


+ H 3 CO+ + H 2 S 










595 


S + H 3 S+ -» 


H 2 S+ + H 










597 


S+ + H 2 S -» 


S+ + H 2 


6.40e- 


10 


-0.50 


0.00 


603 


H 2 S + SO+ - 


-> S+ + H 2 


1.10c- 


09 


0.00 


0.00 


614 


H+ + H 2 - 


> H 2 0+ + H 


7.30c- 


09 


-0.50 


0.00 


624 


H+ + SO -» 


SO+ + H 


1.40c- 


08 


-0.50 


0.00 


739 


S + HS+ -> HS + S+ 










740 


S + H 2 S+ -» 


H 2 S + S+ 










751 


H 2 + CH+ - 


> CH+ + PHOTON 










756 


O + SO -> S0 2 + PHOTON 


3.20e- 


16 


-1.60 


0.00 


757 


CO + S -» OCS + PHOTON 










789 


H 3 0+ + e- 


^ OH + H + H 










798 


H 3 CO+ + e~ 


-> CO + H + H 2 










799 


H 3 CO+ + e~ 


HCO + H + H 










800 


H 3 CO+ + e~ 


-> H 2 CO + H 










803 


CH5O+ + e~ 


-> CH4O + H 










804 


CH5O+ + e~ 


-> H 2 CO + H 2 + H 










812 


H 3 S+ + e~ 


HS + H + H 


1.00c- 


07 


-0.50 


0.00 


818 


H 3 CS+ + e~ 


-> CS + H 2 + H 










819 


H 3 CS+ + e~ 


-» H 2 CS + H 










820 


SO+ + e" -» 


s + 
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825 


HOCS+ 


+ e" 


-> OCS + H 






• 


• < 


» • 



The above analysis of "important" reactions refers to only one perturbation ampli- 
tude. In general, since the equations are not linear, the amplitude of the perturbation 
may influence the results in a non-linear way To check for non-linearity we also ran the 
case where each reaction rate is twice as large as the "standard" one (i.e. a pertubation 
amplitude of 1.0, which doubles the rate of reaction), and still obtained linear variations, 
so that normalized variations are independent of amplitude. One exception concerns CS, 
for which several reactions become important in the latter case. These reactions, listed 
in Table 01 for composition A and Model 2, must be added to the ones shown in Fig. ^ 
for CS. 

Because rate coefficients are often dependent on the temperature, a change in this 
parameter can affect which reactions are important. In particular, an increase of the tem- 
perature to 200 K makes some of the neutral-neutral reactions more important (Table |SJ) 
and some reactions (579, 789, 798, 800, 804 and 825) of Table negligible. On the con- 
trary, an increase of the H2 density to 10 8 cm -3 does not change the results significantly. 

Finally, from Tables 01 and and Fig. ^ we can determine the reactions producing 
the largest variations for the S-bearing species under a set of relevant conditions and a 
reasonable time (10 4 yr) for our models: reactions 90, 153, 190, 191, 740, 756 and 757. 
Reaction 90 is the cosmic ray ionization rate of H2, which is obviously important for 
starting the ion-molecule chemistry. Reactions 153 and 190 are neutral-neutral reactions 
important only at the higher temperature considered (200 K). The other reactions (cf 
Tabled are a collection of ion-molecule (740), neutral-neutral (191), and radiative asso- 
ciation (756 and 757) processes. The importance of these reactions would probably have 
been overlooked had this analysis not been done. It will therefore be crucial to know 
the rate coefficients of the mentioned reactions with high precision. Of these reactions, 
only the rate coefficient for 153 is well determined in the laboratory, although the lowest 
measured temperature (300 K) means that the rate coefficient in the 100-200 K range 
involves an extrapolation. Uncertain measured activation energies for reactions 190 and 
191 also lead to poorly determined rate coefficients by 100-200 K. The ion-molecule re- 
action (740) may not even be exothermic, while the rate coefficients for the radiative 
association processes are order-of- magnitude estimates at best. 
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Table 4. Additional important reactions for a perturbation amplitude of 1.0, for Model 
2 and Composition A. 



Number 


Reaction" 






Q 








7 


168 


C + SO2 — ► CO + SO 
















173 


CH 2 + O^CO + H + H 


1 


20e- 


10 





00 





00 


308 


He+ + H 2 -> OH + He + H+ 


1 


32e- 


09 


-0 


50 





00 


309 


He + + H2O OH+ + He + H 


1 


32e- 


09 


-0 


50 





00 


333 


He+ + OCS -> CS + 0+ + He 


8 


40e- 


10 


-0 


50 





00 


360 


C+ + H 2 HOC+ + H 


1 


80e- 


09 


-0 


50 





00 


444 


CH+ + S -> HCS+ + H 2 
















447 


CH+ + SO -> HOCS+ + H 2 


4 


20e- 


09 


-0 


50 





00 


459 


CH4 + S+ -> H 3 CS+ + H 


1 


40e- 


10 





00 





00 


509 


O + HCS+ S + HCO+ 


5 


OOe- 


10 





00 





00 


510 


O + HCS+ OCS+ + H 


5 


OOe- 


10 





00 





00 


552 


H 3 0+ + H2CO -> H 3 CO+ + H2O 


2 


60e- 


09 


-0 


50 





00 


553 


H 3 0+ + CH4O -> CH 5 0+ + H 2 
















646 


He+ + H 2 -> H 2 0+ + He 


1 


32e- 


09 


-0 


50 





00 


663 


C+ + s -» S+ + c 

















a Only rate coefficients that differ by more than 10 percent from the UMIST value are included. 
Table 5. Additional important reactions at 200 K, for Model 2 and Composition A. 



Number 


Reaction" 




a 




p 


7 


107 


CO + CRPHOT -> O + C 


1 


30e-16 





00 


0.00 


125 


SO + CRPHOT -> SO+ + e" 


1 


30e-14 





00 


0.00 


153 


H 2 + O OH + H 












154 


H 2 + OH -> H 2 + H 


8 


40e-13 





00 


1040 


166 


C + SO S + CO 


3 


50e-ll 





00 





173 


CH 2 + O^CO + H + H 


1 


20e-10 





00 





186 


O + OH -> O2 + H 












190 


O + H 2 S ->■ HS + OH 


9 


22e-12 





00 


1800 


204 


OH + SO -> S0 2 + H 












309 


He + + H2O — ^ OH+ + He + H 


1 


32e-09 


-0 


50 


0.00 


360 


C+ + H 2 HOC+ + H 


1 


80e-09 


-0 


50 


0.00 


436 


CH+ + O HCO+ + H 2 


2 


05e-10 





00 


0.00 


437 


CH+ + O HOC+ + H 2 












646 


He+ + H 2 H 2 0+ + He 


1 


32e-09 


-0 


50 


0.00 


663 


C+ + s S+ + c 













a Only rate coefficients that differ by more than 10 percent from the UMIST value are included. 
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3. Results 

3.1. Abundances with respect to H2 




123456 123456 123456 

Log time (yr) Log time (yr) Log time (yr) 



Fig. 2. Evolution of the SO2, SO, H2S, OCS and CS abundances with respect to H2 as 
a function of time, for the three pre-evaporation gas phase compositions A (left panels), 
B (central panels), and C (right panels) and the four grain mantle compositions (1 to 4, 
from the top to the bottom). The gas temperature is 100 K and the H2 density is 10 7 
cm~ 3 . 

We have run the four models listed in Tabled each with the three different gas-phase 
compositions (A, B, and C) prior to evaporation, for gases with temperatures of 100 K 
and 300 K, and densities of 10 5 cm -3 , 10 6 cm -3 and 10 7 cm -3 . In this section, we give 
some sense of the major features of the chemistry of the sulphur-bearing species. 

Fig. |21 shows the evolution of the abundances of the SO2, SO, H2S, OCS and CS 
molecules, with respect to H2, for the four models and three gas-phase compositions 
(with T=100 K and n(H 2 )=10 7 cm~ 3 ). There are clearly differences among the models, 
at times earlier than 10 4 yr, due to the different initial compositions whereas after 10 4 yr 
the three gas phase compositions give similar results. For example, for Models 1, 2 and 
4, the compositions B and C give lower SO and SO2 abundances than the composition 
A at times earlier than 10 4 yr, whereas the three gas phase compositions give similar 
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Fig. 3. Evolution of the SO2, SO, H2S, OCS and CS abundances with respect to H2 
as a function of time, using pre-evaporated gas phase composition A and the four grain 
mantle mixtures (1 to 4, from the top to the bottom). The gas temperature is 100 K and 
the density is 10 5 cm" 3 (a, left panels), 10 6 cm -3 (b, central panels) and 10 7 cm" 3 (c, 
right panels). 



results for Model 3. For sake of simplicity, in the following, we will discuss composition 
A, but the results do not change qualitatively assuming the B or C compositions. 

Figs.[3]Riid0]show the evolution, at 100 K and 300 K respectively, of the abundances of 
the main S-bearing species for the four grain mantle compositions of Table [3 at the three 
different densities. In all the models depicted, SO and SO2 are the most abundant species 
at late times. The final large amounts of SO2 are more noticeable in Models 2-4, where 
large amounts of gaseous sulphur are available. At a temperature of 100 K and a density 
between 10 6 and 10 7 cm" 3 , the SO2 molecule becomes more abundant than SO after 10 4 - 
10 5 yr, mainly because of the neutral-neutral reaction O + SO — > SO2 + photon. Note 
that this radiative association reaction is critical because of the high abundance of atomic 
O in the pre-evaporative gas. At 300 K, the SO2 molecule is formed less efficiently via this 
mechanism since it possesses an inverse dependence on temperature (see Tabled reaction 
756) but it is still as abundant as SO after 2 x 10 5 yr for a density of 10 5 cm -3 , and after 
10 3 yr for a density of 10 7 cm" 3 (see Fig.|2J|. At 300 K, OH is quickly (< 10 2 yr) formed 
through the reaction H2 + O — > OH + H so that SO2 can be formed by the reaction 
OH + SO — » SO2 + H. Here, the presence of a large atomic oxygen abundance in the 
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1 23456 71 23456 71 234567 
Log time (yr) Log time (yr) Log time (yr) 



Fig. 4. Same as Fig. [3] but with a temperature of 300 K. 



pre-evaporated gas-phase is crucial to produce the hi gh abundance o f OH at early times, 
contrary to what was found in previous models (e. g. Charnlevlll997 h In Model 3, where 
the initial sulphur is mostly in S2, the SO molecule is very quickly (< 10 2 yr) formed, as 
S2 directly leads to SO through the reaction S2 + O — > SO + S. 

Now, let us look at the chemistry of hydrogen sulphide, OCS, and CS. The initial 
H2S abundance (see Table |3J, dips sharply after 10 4 yr in all models but increases after 
10 6 yr for models 2 to 4 at 300 K. The decrease of the H 2 S abundance at 10 4 yr is 
produced by the attack on H 2 S by H 3 + , more abundant than in regions, such as 
hot cores, where water has a large abundance. The reaction between H3 + and H2S yields 
protonated hydrogen sulphide (H3S + ), which dissociatively recombines with electrons to 
form HS at least part of the time: H 3 S+ + e~ -> HS + H + H. The HS product is itself 
depleted by the reaction HS + O -> SO + H. At 300 K, H 2 S is efficiently formed at 
later times (> 10 5 yr) by a series of reactions that starts with the destruction by cosmic 
ray-induced photons of SO and SO + to produce atomic sulphur. Atomic sulphur is then 
hydrogenated into HS first, and then H2S (with intermediate steps in which the unusual 
species HS^, H3SJ, HS2 and H2SJ are formed). Note that at lower temperatures S is 
oxygenated rather than hydrogenated, whereas at 300 K atomic oxygen goes into water 
and, therefore, S can be hydrogenated eventually. 

As in the case of H2S, the initial (adopted) abundance of OCS is that derived from 
observations of this species in the solid state (see Section 2.1 and Table 2). Once in the 
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Table 6. Maximum sensitivity of the abundance ratios to assorted parameters. 



Ratio Changing parameter MRV Set of parameter 



SO2/SO Time 

Temperature 

Density 

Model 

SO2/H2S Time 

Temperature 

Density 

Model 

OCS/H2S Time 

Temperature 

Density 

Model 



16 100 K, 10 7 cm -3 , Mod. 3 

6 10 4 yr, 10 5 cm -3 , Mod. 4 

125 10 4 yr, 100 K, Mod. 3 

20 10 3 yr, 100 K, 10 5 cm" 3 

330 300 K, 10 5 cm -3 , Mod. 2 

630 10 3 yr, 10 7 cm -3 , Mod. 4 

5000 10 3 yr, 300 K, Mod. 4 

3150 10 4 yr, 300 K, 10 7 cm" 3 

6 100 K, 10 7 cm" 3 , Mod. 2 

3 10 4 yr, 10 7 cm" 3 , Mod. 2 

10 10 4 yr, 100 K, Mod. 2 

30 10 4 yr, 100 K, 10 7 cm" 3 



gas phase, OCS is destroyed later than H2S. Under some conditions, e.g. T = 100 K and 
n (H2) = 10 6 -10 7 cm -3 , this molecule maintains a sizable if reduced abundance. Actually, 
at T = 100 K and n(H2) = 10 6 -10 7 cm -3 , there is even a temporary increase in the OCS 
abundance from ~ 10~ 7 to ~ 10" 6 (see Figure|3J). The OCS molecule is mostly formed 
by the radiative association reaction CO + S — > OCS + PHOTON and this reaction is 
aided by the high density. 

Unlike the other species, the CS molecule does not start on grain mantles in our cal- 
culations, but is part of the pre-evaporative gas. For composition A, its initial abundance 
is rather low. At 100 K, the initial CS is destroyed increasingly efficiently as density in- 
creases and its abundance never goes over 10" 9 . At 300 K, CS is efficiently produced at 
high density after 10 3 yr and it can be as abundant as 10~ 8 . For compositions B and C, 
there is significantly more CS present initially, especially for composition B, which con- 
tains 50 times more CS than in composition A. The initial CS is more slowly destroyed 
using compositions B and C than using composition A under the conditions depicted in 
Fig. |21 for Models 2 and 4. 

3.2. Abundance ratios 

Abundance ratios between sulphur-bearing species are exceedingly important because 
the use of fractional abundances (X(i) = N(i)/N(H2)) as chemical clocks introduces an 
additional parameter in the analysis - the H2 column density (N(H2)) in the emitting 
region - which is rarely well constrained observationally. Moreover, the abundance ratios 
are less sensitive to the initial amount of sulphur species compared with absolute abun- 
dances. For example, an initial abundance of S or S2 two times less than those used in 
Model 2 or 3 does not change significantly the abundance ratios. This is why observed 
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Fig. 5. Evolution of the abundance ratios SO2/SO (upper panels), SO2/H2S (middle 
panels) and OCS/H2S (lower panels) as functions of time for composition A, a gas tem- 
perature of 100 K, and densities of 10 5 cm -3 (left panels), 10 6 cm~ 3 (central panels) and 
10 7 cm~ 3 (right panels). The solid, dashed-dotted, dashed, and dotted lines represent 
results from Models 1,2,3, and 4, respectively. In the upper panels, the results of Model 
2 and and 4 are the same and represented by dashed-dotted lines. 



abundance ratios are typically used to put contraints on chemical models, when the two 
species trace the same region. 

Figs. [U and show the evolution of the abundance ratios SO2/SO, SO2/H2S and 
OCS /H 2 S, at 100 K and 300 K for three densities using composition A and all four models. 
The overall sense of these figures is that the computed ratios are more sensitive to the 
gas temperature, the density, and the (poorly known) form of sulphur in grain mantles 
than to the time (particularly for t < 3 x 10 4 yr) or to the gas phase chemical history of 
the cloud. As a first step in quantifying the sensitivity of the studied abundance ratios 
to the different parameters, we report in TableEJthe maximum relative variation (MRV) 
for the three abundance ratios. These are obtained by changing the parameters time, 
temperature, density, model within the following considered ranges: time, 10 3 -10 4 yr; 
temperature, 100-300 K; density, 10 5 -10 7 cm -3 ; grain mantle composition, Models 1-4. 
The MRV is determined for each parameter by varying this parameter within the stated 
range while holding the others fixed at the set of values shown in TableEl These sets lead 
to the strongest changes; hence the term "maximum relative variation". For example, 
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the maximum variation of the SO2/SO ratio produced by an increase of the time from 
10 3 to 10 4 yr is 16 and it is obtained for a temperature of 100 K, a density of 10 7 cm~ 3 
and Model 3. 

The SO2/H2S ratio shows the largest variations with time, but, unfortunately, also 
with differing physical conditions and mantle mixtures. Indeed the ratio is more sensitive 
to the latter parameters than to the time, for t < 3 x 10 4 yr. The OCS/H2S ratio is less 
sensitive to the different mantle mixtures and physical conditions, except at late times 
(> 10 5 yr). Unfortunately, however, this ratio is not sensitive to the time at t < 10 5 yr 
either. The SO2/SO ratio shows significant variations with respect to the different mantle 
compositions, and different densities and temperatures at relatively early times (< 10 4 
yr). These variations mask completely the dependence on the time, which is, anyway, 
moderate at 100 K with the exception of the high density case, and small at 300 K. 

To gain an understanding of the sensitivities of these abundance ratios in addition 
to that obtained from the MRV analysis, we follow the variations of SO2/SO, SO2/H2S 
and OCS/H2S compared with reference ratios computed from Model 2 at 10 3 yr for a 
temperature of 100 K and a density of 10 5 cm~ 3 (left panels of Fig. [SJ. The reference 
SO2/SO ratio is ~ 5x 10 -3 and increases by only a factor of 2 by 10 4 yr. However, a similar 
increase can also be due to an underestimate of the density by a factor of 10 (note that 
if nn 2 = 10 7 cm~ 3 , the SO2/SO ratio increases by a factor of 20). Also, the reference 
SO2/SO ratio at 10 3 yr varies by a factor of ~20 depending on the adopted mantle 
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mixture. The reference SO2/H2S ratio is 0.01 and increases by one order of magnitude 
at 10 4 yr. An increase in density (to 10 7 cm~ 3 ) or in temperature (to 300 K) gives 
respectively a ratio of 0.3 or 0.1. The largest variation is seen if both the temperature 
and the density increase, so that the SO2/H2S ratio jumps to 160. A different mantle 
mixture leads to ratios between 3 x 10~ 3 and 0.3 (at the reference time). The reference 
OCS/H2S ratio is 1 and shows little to no sensitivity to time until 10 4 yr, but also evinces 
little change with temperature, density and mantle composition, with the exception of 
Model 4 which gives a ratio enhanced by a factor of ~10 because the abundance of H2S 
is 10 times less than in Models 1, 2, and 3. 



4. Discussion 



4.1. Comparison with previous models 

The most recent p apers focu s sing o n the detailed mode ll ing of sulphur chemistry in hot 



cores are those of 



Charnlevl 1 1997ft and 



Hatc hell et al 



11998). The adopted chemistry 



is roughly the same in the two models, but they differ in the value of the activation 
energy bar rier for the dest ruction of H2S by atomic hy drogen (Ho S + H — ► HS + H2). In 



particular, 



Hatc hell et al 



used a value of 350 K while 



Charnlevl used a value of 850 K. 



Although there is some experimental support for a variety of choi ces, here we are us ing 



Peng et al 



199fjl) 



an energy barrier of 1350 K, based on the recent experiment by 
the 300-600 K range. The difference in the activation energy barrier for the H2S + H 
reaction causes significant variations in the models. While we obtain results similar to 
the Charnley model for 100 K, H 2 S destruction occurs at later times in our model for 
T = 300 K. Moreover, the high value of the energy barrier has significant consequences 
on the SO and SO2 abundances of our model with respect to the above two models. For 
exa mple at 10 4 yr, our model predicts SO and SO2 abundances 6 and 30 times lower 



than 



Charnlevf s model, using the same parameters and initial conditions. Overall, and in 



addition to this, our results are different from these two models primarily because we are 
assuming different initial gas phase and/or mantle compositions. As already noted (in 
section 3.1), the large atomic oxygen abundance is, for example, an important difference 
in t he pre-evaporat ed gas phase composition. Another difference is that both 



and 



Hatchell et al 



Charnlcv 



assume that the bulk of sulphur in the ice mantle is in the H2S form. 
Our model 1 adopts a similar mantle mixture, with solid-phase OCS as abundant as H2S, 
whereas models 2, 3 and 4 assume that the bulk of the sulphur lies in the S or S2 forms. 
These two classes of models give very different results, as widely discussed in section 3. 

Actually, a novelty of this work is the consideration that the sulphur can be in the 
atomic or molecular form when evaporated from the grain mantles. However, we should 
emphasize that atomic S is very quickly locked into SO and SO2 molecules, so that we 
do not expect large atomic S quantities in the gas phase for a long time. For example, 
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for a gas at 300 K and 10 7 cm -3 , the S abundance decreases fr om 3 x 10~ 5 to 1.5 x 10~ 9 
in 10 3 yr by reacting with O2 to form SO. This explains why Ivan der Tak et alJ (2003) 
did not detect atomic sulphur 2 in the hot core region of the massive protostars that they 
studied, where the gas temperature and density are similar to those quoted above. In a 
slightly colder gas (~ 100 K), atomic S would survive longer, but it would be very difficult 
to detect. In fact the intensity of the SI fine structure line at 25.249 /im is 1 x 10~ 12 erg 
s _1 cm -2 , for a gas at 100 K and 10 7 cm -3 , a source size of 2", an H2 column density 
of 10 23 cm~ 2 , and all the sulphur in atomic form (as our models 2 and 4 adopt). The 
signal would be attenuated by the foreground dust of the cold envelope, whose H2 column 
de nsity is also a round 10 23 cm -2 , by about a factor of ten (using the extinction curves 



Drain putting the signal at the limit of the ISO-SWS detection (a typical rms 



is few times 10 13 erg s 1 cm 2 ). Future investigations in the ISO database will be done 
in order to check the possible presence of atomique sulphur in hot cores. 



Hatc hell et al 



1 1998) found the cosmic ray rate to be an important parameter chang- 
ing the timescales of the destruction of H2S and the formation of SO and SO2. The 
standard value usually assumed is 1.3 x 10~ 17 s _1 but there are some indications that it 

1 : — in 

can vary depending on the region (see Ivan der Tak fc van Dishocck 2000). To check the 



dependence on ionization rate whether or not it can pertain to objects as dense as hot 
cores, we have run the models for three different cosmic ray ionization rates: 1.3 x 10~ 17 , 
1.3 x 10~ 16 and 1.3 x 10~ 15 s -1 . The results for the S-bearing abundance ratios do not 
strongly depend on this parameter at a density of 10 5 cm -3 . At higher densities, how- 
ever, an increase in the cosmic ray ionization rate speeds up the photodestruction of 
OCS, H2S and SO2 via cosmic ray-induced photons. For a cosmic ray ionization rate of 
1.3 x 10~ 15 s , a temperature of 100 K and a density of 10 7 cm~ 3 , the abundances of 
H2S and OCS start to decrease before 10 3 yr and the maximum abund ance of SO2 is 
4.5 x 10~ 8 instead of 2 x 10~ 7 with Model 1 and composition A. Hence, as 



Hatchell et al 



1 1998), we found that the adopted value of the cosmic ray rate may be an important 
parameter at high density whereas it has weak consequences on the results at lower 
densities. 



Recently, iNomura fc Millar! l|2003fl reported a study of the chemical composition 
across the envelope of the massive protostar G24.3+0.15. In this study, they derived 
the density and temperature profile and computed t he chemical compo s ition of the gas 



as a function of the radius and time. Evidently, the 



Nomura fe Millarl l|2003|) model is 



more complete than ours in dealing with the physical structure of the protostar, for this 



van der Tak et al 



I 2003) argued that the non detection of the 25 /im line corresponds to 
an upper limit to the atomic sulphur of < 5 x 10~ 8 , computed by considering line absorption. 
Indeed, this is an upper limit on the S abundance in the absorbing gas, i.e. the foreground cold 
gas, rather than the hot gas. The upper limit on the emitting hot core gas is much higher, as 
explained in the text. 
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is the focus of their model. On the other hand, given its complexity, the model does not 
explore the robustness of the achieved results as a function of the necessary assumptions 
of the mod e l itsel f, which is, on the contrary, the focus of our study. In the same way, 



Dotv et al 



I 2002) studied the chemical composition across the envelope of the massive 



protostar A FGL 2591. Finally, a furthe r degree of complexity has been added to the 



problem bv iRodgers fc Charnlevl {2003ft . who considered the evolution of a protostar, 
including the evolution of the thermal and physical structure of the envelope plus the 
chemical evolution. As in the previous case, the advantage of having a better descrip- 
tion of the protostar comes along with the disadvantage of a lack of exploration of the 
robustness of the results as a function of the chemical network utilized. In these three 
models, the authors assumed th at sulphur was on l y evap orated from the grains mantles 



in the PCS and H2S forms for 



Dotv et al 



Nomura & Millar (2003) and only in the H2S form for 
I20U2) and IRodgers fc Charnle J 1 2003 ) . We have discussed above the prob- 



lem with the assumption that sulphur is totally frozen onto OCS and H2S, and showed 
that the resulting S-bearing abundances depend dramatically on this assumption. The 
conclusion is the same if considering a model where sulphur is only in H2S (as discussed 
at the end of Sect. 4.2). It comes as no surprise, therefore, that our results are different 
from those of these more detailed approaches. 



4.2. Comparison with observations 



In the following, we consider the case of two well studied hot cores: the massive Orion 
KL and the low mass IRAS16293-2422 hot cores. Orion KL is a complex region where 
several energetic phenomena are present. High resolution, interferometric observations 
have shown that different molecules originate in different components, especially for the 
OCS, SO and SO2 molecules (e.g. 



Wright et al 



1996). It is evidently a very crude 



approximation to attribute all the S-molecule emission to the hot core, but lacking a 
better understanding of the region we will try to compare the observed abundances with 
our model predictions. First, in order to minimise the number of parameters, we will 
compare available observations with computed abundance ra t ios. T he S O, SO2, H2S 



and O CS column densities have been derived by 



Schilk e et al 



(2QQl|) and 



Sutton et al 



(1995) to be 2.3 x 10 17 , 6 x 10 16 , 1.2 x 10 16 cm" 2 and 9 x 10 15 cm" 2 , respectively. All 



these column densities are beam-averaged, with a beam of ~ 10". For simplicity, we 
will assume that t he emission origina tes in the hot core region with T = 200 K and 
n(H2) = 10 7 cm~ 3 I Wright et al.lll992 ). The comparison of the abundance ratios derived 



from the above column densities with our model predictions shows that none of the four 
models of Sect. 2.1 reproduces the observed data. However, the ratios can be reproduced 
by a model similar to Model 2, where the initial abundance of atomic sulphur injected 
into the gas phase is 3 x 10~ 6 , i.e. ten times lower than the abundance used in Model 2. 
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Table 7. Observed and modelled abundances of SO, SO2, H2S, OCS and CS towards 
the IRAS16293-2422 hot core. 



Species 


observed 


Model 2 


Model 2' 


SO 


1.7 x 10" 6 


5 x 10~ 7 


4 x 10~ 7 


SO2 


5.4 x 1CT 7 


lCT 7 


10- 7 


H 2 S 


5.3 x 1(T 7 


1CT 7 


10- 7 


OCS 


1 x 1(T 6 


2 x 10~ 7 


lCT 7 


CS 




1 x ltr 11 


2 x 1(T 12 



We will call this model Model 2'. In that case, we obtain a good agreement between the 
model and the observations and derive an age of 4 x 10 3 yr. Assu ming the emission r egion 



1995), the 



extends about 10" and an H2 column density of 8 x 10 23 cm" 2 (Sut ton et al 
observed abundances of SO, SO2, H2S and OCS are reproduced to within a factor of 10. 
On the contrary, the CS abundance is underestimated by three orders of magnitude by 
our model, probably because CS emission does not originate in the hot core. Note that, 
given the complexity of the region, a better analysis would require an understanding of 
the H2 column density associated with the S-bearing molecules. 

In the case of the low-mass protostar IRAS16293-2422, multifrequency observations 
of H2S, OCS, SO and SO2 have been reported and analyzed by sophisticated models, 
which take into account the density and temperat ure structure of the source, as well 



as th e abun dance profile of each stu died molecule dSc hoicr 



2004). Both 



Schoier et al 



(2002) and 



Wakelam et al 



et al 



2002; 



Wak elam et al 



2004) found the following abun- 



dance ratios in the inner warm region where the dust mantles evaporate (jCeccarelli et al 
200o4 : SO 2 /SO=0.3, S0 2 /H 2 S=1 and OCS/H 2 S=1.9. 



The modelling of t he density and temperature profiles of IRAS16293-2422, by 



Ceccarelli et al 



(2000a|), suggests a density of 10 7 cm" 3 and a temperature of 100 K 
in the hot core region. Fig. shows the comparison between the observed and predicted 
ratios (assuming composition A) for the four models: only Model 2 reproduces the obser- 
vations, suggesting an age of ~ 2 x 10 3 yr for the protostar. To check the robustness of 
this result, we also tried a variety of intermediate mantle mixtures, between Model 2 and 
Model 3, with the following characteristics: (1) one-half of the initial mantle sulphur in 
atomic form and one-half in the form of S2, (2) one-quarter in S and three-quarters in S 2 , 
(3) three-quarters in S and one-quarter in S2. In all three intermediate cases, the resulting 
curves are very similar to those found for Model 3. The presence of molecular sulphur 
clearly determines the form of the curves and is inconsistent with these observations. The 
observed ratios are, however, well reproduced by Model 2' (used for the comparison with 
Orion) giving an age, ~3x 10 3 yr, very similar to the age derived by Model 2. Finally, 
both Model 2 and Model 2' predict fractional abundances at the optimum times in good 
agreement (within a factor of 5-10) with those observed, as reported in Tabled Note 
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Fig. 7. Comparison between the observed ratios SO2/SO, SO2/H2S and OCS/H 2 S to- 
wards IRAS 16293-2422 and the theoretical predictions of four models (Mod. 1, 2, 3 and 
4; all with composition A) for a gas temperature of 100 K and a density of 10 7 cm~ 3 . 
The grey empty boxes, if shown, represent those instances in which particular calculated 
ratios are in agreement for a range of times with the observed ones taking into account 
the uncertainties in the observations. The grey filled box for Model 2 shows the time 
interval where all calculated ratios are in agreement with the observed ones. 



that the CS abundance towards IRAS16293-2422 has not been derived 1 Schoier et al 



2002 



Wakelam et al 



2004) and that the abundances of the main S-bearing species pre- 
dicted by Model 2 and 2' are very similar until 10 4 yr. This comparison suggests that 
the majority of the sulphur is released into the gas phase in its atomic form or quickly (t 
< 10 3 yrs) converted to it, and that the abundance of the H2S molecule injected in the 
gas phase from the grains mantles cannot be much less than 10~ 7 . 

The age determined for I RAS16293-2422 is relat i vely short compar e d with the pre 



vious estimates (~ 10 4 yr) l|Ceccarelli et al 



2000b; 



Maret et al 



2002; 



Wakelam et al 



2004). However, our newly determined age represents only the time from the evaporation 



of the mantles; the dynamical time needed to reach the required luminosity to form the 
observed hot core is an additional ~ 3 x 10 4 yr. 

The discussion of the best chemical age for IRAS 16293-2422 is based on models with 
the standard value of the cos mic ray ionization r ate (1.3 x 10~ 17 s~ 4 ). If we use a rate 10 



times higher as suggested by 



Dotv et, al 



the observations towards IRAS 16293- 



2422 are no longer in agreement with Model 2, but with Model 1, although the derived age 
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is very similar. Indeed, the enhanced cosmic ray ionization rate speeds up the destruction 
of H2S, increasing the ratio SO2/H2S more rapidly without affecting significantly the 
other ones. Consequently, the SO2/H2S curves on Fig. [7| are shifted to the right of the 
figures, worsening the agreement with Model 2 but improving it with Model 1. In that 
case, the derived age is 5 x 10 3 yrs and the predicted absolute abundances are between 15 
and 20 less than the observed ones, however. Thus, a high cosmic ray ionization rate is 
no longer compatible with our hypothesis that mantle evaporation leads quickly to high 
gas-phase abundances of atomic sulphur. The high cosmic ray ionization rate has also a 
consequence on the O2 abundance pedicted to be ~ 10~ 6 at 10 3 yr with Model 1 whereas 
it is predicted to be two orders of magnitude lower with the standard value of the cosmic 
ray ionization rate and Model 2. However, at present one cannot choose between the two 
conditions (large abundances of S in the gas phase or large cosmic-ray ionization rate). 
High sensitivity observations of atomic sulphur (see Sect. 4.1) are needed to put stringent 
constraints on chemical models. 



Another way to confirm our hypothesis of a large initial abundance of atomic sulfur 
would be a careful study involving several other protostars at different stages of evolution. 
The strongest prediction of our models 2 to 4 is the large (<; 10 -5 ) amount of SO2 at late 
times, compared with models that start with little sulphur evaporated from the grains. 
There is some evidence for SO2 abundances in high-mass hot cores as large as 1CP 6 



I van der Tak et al 



2003), suggesting that an initial amount of sulphur at least higher 



than 10~ 6 is needed. But our prediction of large SO2 abundances in more evolved hot 
cores with ages larger than ~ 10 5 yr, does not pertain to the Orion and IRAS16293-2422 
hot cores (as found in the present work). 

Quantitative comparisons with o bservations of other hot c ores reported in the lit — 



crature, such as the observations by 



Hatchcll et al 



(1998) and 



Buckle fc Fulled (|2003), 



are difficult to carry out since the abundances determined by these authors are averages 
along the line of sight and beam-averaged. He nce, their analyses do not take into account 



van der Tak et al 



1 2003) report the study 



the physical structure of the source. Finally, 
of the S-bearing molecular abundance in half a dozen massive protostars. In this case, an 
attempt to disentangle the outer envelope and inner core has been done, but the analysis 
is not accurate enough to derive abundance ratios to compare with detailed model predic- 
tions; rather, th e authors just give order- of-magnitude jumps of abundance in the inner 



hot core regions. 



van der Tak et al 



model predictions by 



Dotv et al 



2003) compared the measured abundances with the 



2002), and argued that OCS is the main sulphur bear- 
ing molecules on the mantles, because its predicted abundance is otherwise too low. We 
ran a model where sulphur is in solid OCS, to test this suggestion, and we found that it 
leads to results very similar to those of Model 1 (where sulphur is evaporated from the 
grain mantles in the H 2 S and OCS forms), because OCS has almost the same chemical 
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behavior as H2S. Although OCS is destroyed at later times, the differences in the SO 
and SO2 formation are not significant and the same discussion about Model 1 is valid. 



5. Conclusions 

We have studied in detail the influence of the mantle form of sulphur on the post- 
evaporative gas-phase abundances of S-bearing molecules in hot star- forming regions, 
with the goal of understanding whether those molecules can be used to estimate ages. 
We considered four different reasonable mantle mixtures, from which gas-phase H2S, 
OCS, S and S2 emerge after a process of evaporation and, for the last two species, 
possible rapid reaction, with different relative abundances, joining other species in the 
gas-phase prior to evaporation. We then followed the post-evaporative chemical evolution, 
with an emphasis on the abundance ratios of the main sulphur-bearing species for realistic 
physical conditions present in hot cores. Our results show that none of the ratios involving 
the four most abundant S-bearing molecules, namely H2S, OCS, SO and SO2, can be 
easily used by itself for estimating the age, because the ratios depend at least as strongly 
on the physical conditions and on the adopted grain mantle composition as on the time. 
Also, the abundance of atomic oxygen in the gas phase, if not correctly accounted for, 
can seriously affect the chemistry. The situation, however, is not totally hopeless, because 
a careful comparison between observations and model predictions can give some useful 
hints on time estimates, and on the mantle composition. Such a careful analysis has 
to be done on each single source, however, for both the physical conditions and mantle 
composition can vary from source to source, so that the abundance ratios are not directly 
comparable. In practice, a careful derivation of the molecular abundances (which takes 
into account the source structure) coupled with a careful modeling of the chemistry at 
the right gas temperature and density is necessary. 

We applied our model to two well studied hot cores: Orion KL and IRAS16293. For the 
S-bearing abundances towards Orion KL, we assumed that their emission arises from the 
hot core region (which is strongly debatable) and is not beam-diluted. We were not able 
to reproduce all of the observed abundances ratios with any of our models. The agreement 
with Model 2 is satisfactory if we decrease the initial amout of atomic sulphur by a factor 
of 10. In that case, we derive a best age of 4 x 10 3 yr. However, the predicted abundance 
of CS is three orders of magnitude lower than the observed one. Contrary to the case of 
Orion KL, the sulphur-bearing abundances though the low mass ho t core of IR AS 16293- 



2422 have been carefully determined through a sophisticated model l(Schoier et allLoO^ . 
which takes into account the density and temperature structure of the source, as well 
as the abundance profile of each studied molecule. Using the standard value of cosmic 
ray rate, we found that Model 2, in which a large amount of atomic sulphur is initially 
present in the post-evaporative gas, best reproduces the observed abundance ratios. In 
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that case, we derived an age of ~ 2 x 10 3 yr from the evaporation era to the current stage 
of this particular low mass hot core. If we decrease the initial amount of atomic sulphur 
in Model 2 as for Orion KL, the agreement is still good and gives a similar age. This 
analysis favors the hypothesis that sulphur is mainly evaporated from the grains in the 
atomic form or in a form quickly convert ed into it. On the c ontrary, if a higher rate is used 



as suggested by the recent modelling of 



Dotv et al 



I 2004J), best agreement occurs with 



Model 1, where no atomic sulphur can be found in the grain mantle and only H2S and 
OCS are initially present. The strongest prediction of our atomic sulfur-rich model is the 
presence of large abundances of SO2 , derived from this form of sulfur, at late stages of hot 
cores. A futher systematic study of S-bearing-species towards older hot cores where the 
physical structure is well known would provide information to test this model. Moreover, 
the fact that not all of the sulfur need be initially in atomic form, given the reasonable 
agreement obtained using Model 2', suggests that a signficant port ion of the granula r 



elemental sulphur may be tied up in materials such as iron sulphide l|Keller et al 



2002). 
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